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Nailing down the unknown neutrino mixing angle #13 is one of the most important goals in current 
lepton physics. In this context, we perform a global analysis of neutrino oscillation data, focusing on 
813, and including recent results [Neutrino 2008, Proceedings of the XXIII International Conference 
qq j on Neutrino Physics and Astrophysics, Christchurch, New Zealand, 2008 (unpublished)]. We discuss 

, two converging hints of 813 > 0, each at the level of ~lcr: an older one coming from atmospheric 

neutrino data, and a newer one coming from the combination of solar and long-baseline reactor 
£\j ' neutrino data. Their combination provides the global estimate 

sin 2 0i3 = 0.016 ±0.010 (Iff), 

implying a preference for #13 > with non-negligible statistical significance (~90% C.L.). We 
discuss possible refinements of the experimental data analyses, which might sharpen such intriguing 
indication. 
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^ (— 1 PACS numbers: 14.60.Pq, 26.65.+t, 28.50.Hw, 95.55.Vj 

Introduction. — In the last decade, it has been established that the neutrino states (v ei v^,u T ) with definite flavor 
are quantum superpositions of states (yi, v-i, V3) with definite masses (mi, m-2, WI3) Q- These findings point towards 
new physics in the lepton sector, probably originating at very high mass scales 

Independently of the origin of neutrino masses and mixing, oscillation data can be accommodated in a simple 
theoretical framework (adopted hereafter), where flavor and mass states are connected by a unitary mixing matrix 
U, parametrized in terms of three mixing angles (#12, #13, #23) and one CP-violating phase 5 []]]. The mass spectrum 
' gaps can be parametrized in terms of Sm 2 = to 2 — m i and of Am 2 = m 2 — (m 2 + m|)/2 Q. 
\Q . Within this framework, the mass-mixing oscillation parameters (5m 2 , sin 2 612) and (Am 2 , sin 2 #23) are rather well 
determined [3j . Conversely, only upper bounds could be placed so far on sin 2 #13 , a dominant role being played by 
. the null results of the short-baseline CHOOZ reactor experiment [3] (sin 2 $13 < few%). 

Determining a lower bound for #13 (unless #13 = for some unknown reason) is widely recognized as a step of 
paramount importance in experimental and theoretical neutrino physics [HQ. Indeed, any future investigation of 
leptonic CP violation (i.e., of 5), and of the neutrino mass spectrum hierarchy [i.e., of sgn(Am 2 )] crucially depends 
on finding a nonzero value for #13. A worldwide program of direct #13 measurements with reactor and accelerator 
neutrinos is in progress, as recently reviewed, e.g., at the recent Neutrino 2008 Conference 

In this context, any indirect indication in favor of #13 > becomes highly valuable as a target for direct searches. 
We report here two indirect, independent hints of #13 > 0, one coming from older atmospheric neutrino data, and 
one from the combination of recent solar and long-baseline reactor data, as obtained by a global analysis of world 
oscillation searches. For the first time, these hints add up to an overall indication in favor of #13 > at non-negligible 
confidence level of ^90%. 

Hint from atmospheric neutrino data. — In a previous analysis of world neutrino oscillation data Q , we found a 
weak hint in favor of #13 > 0, at the level of ~0.9cr, coming from atmospheric neutrino data combined with accelerator 
and CHOOZ data (see Figs. 26 and 27 in [3(). We traced its origin in subleading Zv oscillation terms driven by Sm 2 
61, which are most effective at cos 5 — -1 (see Fig. 24 in [3]), and which could partly explain the observed excess 
of sub-GeV atmospheric electron-like events Such hint has persisted after combination with further long-baseline 
(LBL) accelerator neutrino data @, @|, which have not yet placed strong constraints to v e appearance. In particular, 
after including the Main Injector Neutrino Oscillation Search (MINOS) data [l(| presented at Neutrino 2008 11 1, 
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and marginalizing over the leading mass-mixing parameters (Am 2 , sin #23)1 we still find a ^0.9er hint of #13 > from 
the current combination of atmospheric, LBL accelerator, and CHOOZ data, 

sin 2 6»i3 = 0.012 ± 0.013 (la, Atm + LBL + CHOOZ), (1) 

where the error scales almost linearly up to ^3er, within the physical range sin 2 #13 > 0. 
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FIG. 1: Allowed regions in the plane (sin 2 612, sin 2 #13): contours at la (dotted) and 2a (solid). Left and middle panels: solar 
(S) and KamLAND (K) data, both separately (left) and in combination (middle). In the left panel, the S contours are obtained 
by marginalizing the Sm 2 parameter as constrained by KamLAND. Right panel: All data. 



Hint from solar and KamLAND data. — In past years, the above "atmospheric v hint" was not supported by 
independent long-baseline reactor and solar neutrino data, which systematically preferred $13 = as best fit, both 
separately and in combination Therefore, in the global data analysis, the hint of #13 > was diluted well below 
lcr, and could be conservatively ignored [3j. 

Such trend has recently changed, however, after the latest data release from the Kamioka Liquid Scintillator Anti- 
Neutrino Detector (KamLAND) [12j, which favors slightly higher values of sin 2 #12, as compared to solar neutrino 



data |l3j at fixed 6*13 = 0. As discussed in [14(, and soon after in |l5j], this small difference in sin 2 #12 can be reduced 
for 0i3 > 0, due to the different dependence of the survival probability P ee = P[y e — » v e ) on the parameters (#12, #13) 
for solar and KamLAND neutrinos [la]. Indeed, recent combinations of solar and KamLAND data prefer 6*13 > 0, 
although weakly 0, H3] . 

Remarkably, the recent data from the third and latest phase of the Sudbury Neutrino Observatory (SNO) [l8| pre- 
sented at Neutrino 2008 |19j further reduce the solar neutrino range for sin 2 #12 and, in combination with KamLAND 
data, are thus expected to strengthen such independent hint in favor of #13 > 0. We include SNO-III data in the form 
of two new integral determinations of the charged-current (CC) and neutral current (NC) event rates [18j . with error 
correlation p ~ —0.15 inferred from the quoted CC/NC ratio error [la ], but neglecting possible (so far unpublished) 
correlations with previous SNO data [3] • W e ig nore the SNO-III elastic scattering (ES) event rate [2^ , which appears 



to be affected by statistical fluctuations [IS 



ES rate measured by Super-Kamiokande 2l| . 



19l | and which is, in any case, much less accurate than the solar neutrino 



In the solar neutrino analysis, we u pda te the total Gallium rate (66.8 ± 3.5 SNU) [22| to account for a recent 
reevaluation of the GALLEX data [13, HI]. The latest Borexino data [H, 0], presented at Neutrino 2008 [13], are 



also included for the sake of completeness. We do not include the Super-Kamiokande phase-II results [28j, which 
would not provide significant additional constraints. Finally, concerning KamLAND, we analyze the full spectrum 
reported in [T3 ] . and marginalize away the low-energy geoneutrino fluxes from U and Th decay in the fit. We have 
checked that our results agree well with the published ones (in the case #13 = 0) both on the oscillation parameters 
{8m 2 , sin 2 6>i 2 ) and on the estimated geo-^ fluxes [5§]. 

Figure 1 (left panel) shows the regions separately allowed at lcr (A% 2 = 1, dotted) and 2cr (A^ 2 = 4, solid) 
from the analysis of solar (S) and KamLAND (K) neutrino data, in the plane spanned by the mixing parameters 
(sin 2 612, sin 2 0i3 ). The 8m? parameter is always marginalized away in the KamLAND preferred region (which is 
equivalent, in practice, to set 8m? at its best-fit value 7.67 x 10~ 5 cV 2 ). The mixing parameters are positively and 
negatively correlated in the solar and KamLAND regions, respectively, as a result of different functional forms for 
P ee (sin 2 9i2, sin 2 #13 ) in the two cases. The S and K allowed regions, which do not overlap at lcr for sin 2 #13 = 0, 
merge for sin 2 #13 ~ few x 10~ 2 . The best fit (dot) and error ellipses (in black) for the solar+KamLAND combination 
are shown in the middle panel of Fig. 1. A hint of #13 > emerges at ~1.2cr level, 

sin 2 6>i 3 = 0.021 ± 0.017 (lcr, solar + KamLAND) , (2) 



with errors scaling linearly, to a good approximation, up to ~3cr. 
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FIG. 2: Global v oscillation analysis: Allowed la ranges of sin #13 from different input data. 



Combination. We have found two independent hints of #13 > 0, each at a level of ~1(T, and with mutually consistent 
ranges for sin 2 #13. Their combination reinforces the overall preference for $13 > 0, which emerges at the level of ~1.6cr 
in our global analysis. In particular, Fig. 1 (right panel) shows the la and 2a error ellipses in the (sin 2 612, sin 2 $13) 
plane from the fit to all data, which summarizes our current knowledge of electron neutrino mixing [30J . Marginalizing 
the sin 2 #12 parameter we get 



sin 2 6>i3 = 0.016 ± 0.010 (la, all oscillation data) , 



(3) 



with linearly scaling errors. This is the most important result of our work. Allowed ranges for other oscillation 
parameters are reported separately [3l| . Summarizing, we find an overall preference for 0i3 > at ~ 1.6a or, 
equivalently, at ~90% C.L., from a global analysis of neutrino oscillation data, as available after the Neutrino 2008 
Conference. The preferred la ranges are summarized in Eqs. (1)— (3), and are graphically displayed in Fig. 2. 

Conclusions and Prospects. — In this Letter, we have focused on the last unknown neutrino mixing angle 0i3. Within 
a global analysis of world neutrino oscillation data, we have discussed two hints in favor of 0i3 > 0, each at the level of 
~1<7. Their combination provides an overall indication for #13 > at a non-negligible 90% confidence level. To some 
extent, the present hints of 0i3 > can be corroborated by more refined analyses. Concerning atmospheric neutrinos, 
an official, complete ?>v analysis by the Super-Kamiokande collaboration, including all experimental details, would 
be very important. The analysis should include <5m 2 -driven terms in the oscillation probability [32I 33 1, which have 
been neglected in the official publication 34|. Concerning solar neutrinos, a detailed, fully documented and official 
combination of all the SNO-I, II, and III data [35[ would be helpful to sharpen the bounds on solar v e mixing and to 
contrast them with (future) KamLAND data. The latter would benefit by a further reduction of the normalization 
error, which is directly transferred to the mixing parameters. In our opinion, such improvements might corroborate 
the statistical significance of the previous hints by another ^la but, of course, could not replace direct experimental 
searches for 0i3 at reactors or accelerators. Two hints make for a stronger indication, but do not make for a compelling 
proof. 
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